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1. Introduction 

Orthogonal polynomial ensembles are given by joint probability densities of the form 

n 

p{Xi,...,Xn) = CnYli^i - ^jfYi^i^i) (1-1) 

i<j i=l 

with a positive weight function w. The classical examples, where w is the weight of Hermite, 
Laguerre, or Jacobi polynomials, naturally arise in random matrix theory: as joint eigenvalue 
distribution of the Gaussian Unitary Ensemble (Hermite), of Wishart matrices (Laguerre), or 
of random projectors (Jacobi). 

In [5], Borodin studied more general joint probability densities of the form 

n 

p{xi, ...,Xn) = CnY\_{^i - - a;^) n^i^(a;i), (1-2) 

i<j i=l 

where 6' is a fixed positive number. As pointed out in [5j, these are related to biorthogonal 
polynomials, and will be referred to as biorthogonal ensembles. 

This note is mainly motivated by a special biorthogonal ensemble that arises from a random 
matrix model for disordered bosons that was proposed by Lueck, Sommers, and Zirnbauer in 
|13] . The model amounts to the product of a Wishart matrix and the fundamental matrix 
of the standard symplectic form, and on the level of eigenvalues, interpreted as characteristic 
frequencies of disordered quasi-particles, one obtains the joint density 

n 

Qn,a{Xl, ■ ■ ■ , Xn) = Cn'^^X°'e ™* | Xj — Xj | | — x|| 1 [0,00)" (^^l , • . . , X„) (1-3) 

1=1 i<*<i<»i 

with a G N. In [13] is was shown that the correlation functions of the frequencies in the bulk 
of the spectrum are in the Gaussian Unitary Ensemble universality class, yet a novel scaling 
behaviour is found at the low frequency end of the spectrum. Other applications of biorthogonal 
ensembles to physics are discussed in [2] and [T7]. In the latter reference, the motivation comes 
from matrix models for Chern-Simons theory. 

The aim of this note is to complement the Lueck-Sommers-Zirnbauer results by a large 
deviations principle for the empirical measure of the characteristic frequencies. For any sequence 
of random numbers Xi, . . . ,Xn we denote by 




i=l 



the empirical distribution or empirical measure of these values (a random probability measure 
on M). Define the mean empirical measure L„ = EL„ by the relation (Z^, /) = E(L„, /) for all 
continuous and bounded functions / : R — t- M. One result in [TB] is that the sequence (-L„)„ 
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converges weakly to a probability measure on M with Lebesgue density 

for < t < 6 := 3^3. One consequence of the large deviations principle that will be proven in 
what follows is that this statement can be improved: we will show that the empirical measures 
{Ln)n themselves converge weakly, in probability, to Qoo- 

Actually, we will study large deviations principles for empirical measures in a broader frame- 
work that encompasses not only ensembles like f 1 1.2 1) , but also takes care of weight functions 
Wn, depending on n, and determinantal parts like Y[i<i I 

for any /3 > 0. We 

will also obtain large deviations results for multiple orthogonal polynomial ensembles (see [H]). 
The main observation is that while these large deviations results do depend on the explicit 
formulae for the joint distributions of the eigenvalues like (11. ip and (11 .2^ . they are independent 
of the determinantal structure of the correlation functions, and we need not invoke the theory 
of orthogonal or biorthogonal or multiple orthogonal polynomials. 

In [2] (see also jjL]|), Ben Arous and Guionnet have obtained a large deviations principle (LDP) 
for the empirical measure of eigenvalues from the Wigner-Dyson ensembles GOE/GUE/GSE 
in the space A^i(]R) of probability measures on the Borel sets of R, endowed with the weak 
topology, with speed n"^ and good rate function (GOE case) 

whose unique minimiser is the semicircle distribution. Recall that a family of probability 
measures {^e)e>o on a topological space X is said to obey a large deviations principle (LDP) 
with speed and good rate function / : X — )■ [0, oo] if J is lower semi- continuous and has 
compact level sets Nl := {x E X : I{x) < L}, for every L G [0, oo), and 



for every open G C X and 



liminf elog/i£(G) > — inf I{x) 



limsup£:log/ie(y4) < — inf I{x) 

e->0 ^eA 



for every closed A C X. 

This has been generalized in [8] to joint densities of the form 

g„(xi, . . . ,a;p(„)) = — JJw„(xi)" JJ ' (l-^) 

" «=i i<«<i<p(") 

with 6 E N, P > 0, partition function Z„ and continuous weight functions : IR — ?■ ■ This 
framework takes care of the needs of mesoscopic physics in that it subsumes matrix versions 
of all classical symmetric spaces (see [10]). With regard to the symmetry classification of 
disordered fermionic systems provided in the last reference, let us remark that an analogous 
classification for the case of bosons is not completely understood, see, however, the discussion 
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in fl8\ Section 4]. Moreover, a d-dimensional generalisation of the random matrix model for 
disordered bosons in [13j was studied recently in P^ . 

The note is organised as follows. Section 2 is devoted to the formulation of a large deviations 
principle for the empirical measures (Iv„)„ of the eigenvalues of generalised biorthogonal ma- 
trix ensembles. The examples include the random matrix model of disordered bosons in [13], 
the Stieltjes-Wigert ensembles in [17] as well as biorthogonal Jacobi-, Laguerre- and Hermite 
ensembles considered in |5]. In Section 3 we formulate large deviations principles for a special 
multiple orthogonal ensemble, the Angelesco ensemble, see [TI] . In Sections 4 and 5 we present 
the proofs of our large deviations principles. 

2. Large deviations for biorthogonal ensembles and beyond 

In this section, we will derive a LDP for the bosonic ensemble, where the density of the joint 
distribution of the eigenvalues is of form (11.31) . Obviously, (11. 3p is a special case of the density 

p(") 

QriiXl, . . . , Xp{ri)) = Wn(Xi) l^^j — | | Xj — | l2p(n) (Xi , . . . , Xp(„) ) , (2.1) 

" i=l l<«<i<p{") 

with 6* G N, partition function Z„ and continuous weight functions w„ : M — )■ . For 9 even, 
S is a closed subset of [0, oo) while for 6 odd, E is a closed subset of R. The sequence (p(n))„ 
must satisfy 

vin) 

lim = K G (0,cx)). 

n— >-oo n 

Note that for p{n) = n and weight functions w independent of n, (12. ip subsumes the density for 
the eigenvalue distribution for biorthogonal ensembles as introduced in [S]. For ^ = 1, p{n) = n 

1 2 

and w(x) = e~2^' we also recover the classical GUE. 

Throughout the whole section, we write M{f) for the set of zeros of a function / : M — )■ M 
and we assume that the sequence of weight functions (w„)n satisfies the following: 

(al) there exists a continuous function : S — )■ [0, oo) such that 

— #A/'(u;) < oo, Uiwn) C Miw) for large n. (al.l) 

— As n — )■ oo, w„ converges to w, and \ogWn to \ogw uniformly on compact sets. 
(al.2) 

— \ogw is Lipschitz on compact sets away from Miw). (al.3) 
(a2) If S is unbounded, then there exists rio G N such that 

lim \x\'^S+l){K+e) ^^^^^ ^ Q 

for some fixed e > 0. 
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We will study the asymptotic behaviour of the empirical distribution of x for x = 

(xi, . . . , Xp(„)) G T,P^"''>, which is defined as 

^ p(") 
Ln(x) := , , > 5x , 

Let the space of probability measures on the Borel sets of E be denoted by A^i(S) and let Qn 
denote the joint distribution of random variables (Xi, . . . ,Xp(„)) with density g„. The main 
theorem now reads as follows: 

Theorem 2.1. The sequence {Qn o L~^)^ satisfies a LDP on A^i(S) with respect to the weak 
topology with speed and good rate function 

/(/i) = y/ j {\og\x' -y'\-' + \og\x-y\-'] ^l{dx)^^{dy) (2.2) 
—K j \ogw{x) ^{dx) + 

where /i G A^i(S) and 
c:=lim^logZ„ = - ird f f {\og\x^ - y^\ ^ + \og\x - y\-^\ ^{dx) ^i{dy) 



-n J log w (x) fi{dx)^ < oo. 



(2.3) 



Corollary 2.2. Whenever the joint density of the eigenvalues is as in (12.11) and the rate func- 
tion I has a unique minimiser /x*, we obtain under (al) and (a2) a strong law of large numbers, 

P(L„ ^ /i*) = 1. 

Moreover for any closed A C A^i(S) with Afl {z/ G A^i(S) : I{v) = 0} = 0, then for all n > hq 

Qn{Ln eA)< exp(-n2 inf /(z/)/2). 

Proof. Employing the upper bound of the LDP, the strong law of large numbers follows via an 
application of Borel- Cantelli's lemma, see |9j, Thm. II. 6. 3]. □ 

Example 2.3 (Disordered bosons). Returning to the bosonic ensemble with density (IL3p . we 
have as weight functions 



Wnix) = X"e " . 



Now r^^ is the variance of the independent and normally distributed random variables, that 
were used to construct the stability matrix h for that ensemble, cf. [13]. We will take 
the variance equal to n~^. For a greater generality we take a sequence {a{n))nm with 
lim„_>oo = a > —1. The case a = incorporates having a constant sequence {a{n))n- Ob- 
viously, conditions (al) and (a2) are met and we obtain for the ensemble of disordered bosons. 
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that [Qn o L^^)n satisfies a LDP on A^i(S) with respect to the weak topology with speed ra^ 
and good rate function 

Hf^) = -y J J log {\^^ -y'^\\x- y\) /i(rfx) fi{dy) 
—K J a logx fx{dx) + K J X fj,{dx) + c, 

where /i G A^i(S) and 

c := hm — log Z„ < oo. 

n^oo n 

In [13], the authors proved the weak law of large numbers (WLLN) 



E 



1 ^ 1 /"^ 

]^X^1{A.<^-} / Poo{t)dt ioT N ^ OO, 

i=l J -^""^ 



with ^oo defined in f ll.4p . They also stated that the corresponding rate function / is convex 
and thus, one can find a unique minimiser fi* of J, /(/U*) = 0, with density poo- Therefore the 
WLLN of [13] is extended into a strong one. 



Example 2.4 (Stieltjes-Wigert ensembles). In the case of Stieltjes-Wigert ensembles form [T7] 
the weight function is w{x) = e~'^'^^°^*^^^)^ with some constant c. Hence our Theorem applies and 
we obtain a LDP and a strong law of large numbers. 

Example 2.5 (Jacobi-, Laguerre- and Hermite biorthogonal ensembles). In [5], Jacobi ensem- 
bles with w{x) = on (0, 1) are studied. We consider the general case 

S = [0,1], w;„(x) =x"/"(l-x)^/", a,(3>-l. 

For 6 = 1 these ensembles appear in the canonical correlation analysis, where the correlations 
coefficients turn out to be the square root of the eigenvalues of a special kind of matrix. It is 
known that the eigenvalues of these matrices follow a joint distribution of Jacobi type. Taking 
a, /3 > —1 fixed and — )■ n, the sequence {Qn o L~^)n obeys a LDP with speed and rate 
function (up to a constant) 



Y / Ji^^s\^^~y^\ ^ + '^og\x-y\ ^)ii{dx)^{dy). 



(2.4) 



For — )■ a > — 1 and — ?• /3 > —1, the corresponding LDP holds true, and the rate 
function is (up to a constant) (12. 4p plus 



-K J a log xfi{dx) — K J /31og(l 



x)jj,{dx). 



Hermite ensembles can be considered with Wn{x) = |x|°/"e~^' on (—00,00) with a > — 1. 
The case a = is the classical Hermite weight. For constant a > — 1 we obtain the LDP, and 
the rate function is (12. 4p plus k J x'^ji^dx) with the same — )■ k. For — )■ a > — 1, the rate 
is (12. 4p plus K J x'^fi{dx) — n J alogxfi{dx). Finally, the Laguerre biorthogonal ensembles are 
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given by the weight function = x°'f'^e~^^^ on (0, oo) with a > —1. It is straightforward 

to find the corresponding rate functions for constant and n-dependent parameters. For 6 = 2 
the Laguerre case is corresponding to the matrix model of disordered bosons in |13] . 

Remark 2.6. It is obvious how to generalise biorthogonal ensembles. Consider an ensembles 
of n points on (a, 6) C M with the joint probability density of the form 

n 

cHwnix.r det[Ux,)]l,=, det[r/,(x,)]i:,=i, (2.5) 

i=l 

where ^i{x) , rii{x) , i > I, are some functions defined on (a, 6). In case p{n) = n the density 
( 12. ip is clearly a special case of (12. 5p . take ^j(x) = and rji{x) = x^^''^^\ It is known, 
see [7] and references therein, that random matrix theory provides many instances of such 
biorthogonal structures, for example unitary invariant matrix ensembles or unitary ensembles 
with an external source, see also [B] and the next section. But large deviations principles for 
the corresponding empirical measures of n points distributed according to (12. 5p are out of range 
with respect to the techniques we are using to prove Theorem 12.11 



3. Large deviations for multiple orthogonal ensembles 



Multiple orthogonal polynomials are a generalisation of orthogonal polynomials in which the 
orthogonality is distributed among a number of orthogonality weights. They appear in random 
matrix theory in the form of special determinantal point processes that are called multiple 
orthogonal polynomial (MOP) ensembles. In pTl [12] the appearance of MOP in a variety of 
random matrix models and models related with particles following non-intersecting paths have 
been considered. To a finite number of weight functions Wi,...,Wp on M and a multi-index 
n = {rii, . . . , Up) G W we associate a monic polynomial of degree n := \n\ := rii + ■ ■ ■ + Up 
such that 

poo 

Pfi{x)x^Wj{x) dx = 0, for = 0, . . . , Uj — 1, j = 1, 



If Pa uniquely exists then it is called the multiple orthogonal polynomial (MOP) associated with 
the weights wi, . . . ,Wp and multi- index n. In [IT] the following result was presented. Assume 
that 



det[fj{xk)]j 



fc=i, 



(3.1) 



is a probability density function on . 
linear span of {x''Wj{x)\j = 0, . . . ,nj 



l<j<k<n ^ 

, where the linear span of /i ,...,/„ is the same as the 
1, j = 1, . . . ,p }. Then the MOP exists and is given by 



Pnix) 



E 



Y[{X - Xj) 

Lj=l 

where the expectation is taken with respect to the p.d.f f l3.ip . which can be interpreted as the 
expectation of the random polynomial YYj=ii^ ~ ^j) with roots xi, . . . ,Xn from a determinantal 



8 



PETER EICHELSBACHER, JENS SOMMERAUER AND MICHAEL STOLZ 



point process on the real line. The p.d.f (13. ip is called a MOP ensemble. It was first observed 
in [5] that random matrix models with an external source lead naturally to MOP ensembles. 

The weights wi, . . . ,Wp are an Angelesco system if there are disjoint intervals Fi, . . . , Fp C M, 
such that supp(wj) C Tj, j = l,...,p. In the Angelesco case, det[fj{xk)]j,k=i,...,n is of block 
form and results in 

det[/,(x.)],,,=i,...,. = WU{X^^)Wwlxf) 

i=\ ^ k=\ 

with := a;Ar._^+fc G Fj, A^^ = X]j=i "^j (with A'o = 0) and X^*) = [x"^ , . . . , xi*]) and 

A(X) = JJ {xk-Xj) for X = (xi, . . . ,a;„), 

the Vandermonde determinant. Thus an Angelesco system gives rise to a MOP ensemble, the 
Angelesco ensemble^ and the joint p.d.f is 

i-nA(x«)^ n A(x«,x(^-))nnu;.(x«), (3.2) 

i=l i<*<i<p *=1 fc=l 

where 

A(x,y):=nn(^'^-^^) 

fc=i i=i 

for X = (xi, . . . , x„) and Y = {yi, ...,?/„). 

We now consider the situation that |ri| = n — )■ oo and nj — )■ cxd for every j = 1, . . . ,p in such 
a way that 

— ^ Tj forj = l,...,p (3.3) 

n 

with < r j < 1 and ^j^iT'j = 1- Let us consider varying weights 

Wi{x) = e-"^'(^) (3.4) 

for any i = 1, . . . ,p. Denote by 

Uj 

L,(x(^)):=L,,„^,(x(^)):=-X;5^0) (3-5) 

J k=l 

the j-th empirical measure of x'^^'^ for every j = 1, . . . 

We will study the asymptotic behaviour of the empirical distribution vector Ln{x) of x for 
X = (xi, . . . , x„) G 11^=1 =■ Tp.n, which is defined as 

L„(x) :=(Li(x«),...,Lp(x(P))) 

with Lj(x(*)) defined as in ( 133|) . L„(x) is an element in xf^^A^i(F"*) =: A^i(Fp^^). Let 
denote the joint distribution of random variables (X^^-*, . . . ,X*-^-') with density (13. 2p . 
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Theorem 3.1 (LDP for Angelesco ensembles). Assume that every weight function Wi in (13. 4p 
satisfies assumption (al) and (a2) and assume that (13. 3 p is fulfilled. Then the sequence 
{Qn o L~^)^ satisfies a LDP on M.i{Tp^fi) with respect to the weak topology with speed and 
good rate function 



/(/ii,...,/ip) = ^^r^ I I ^og\x - y\ ^ fij{dx) fij{dy) (3.6) 



p-i p 



=i+i 

where (/^i, . . . , /ip) G AiiiTp^fi) and 



+ XI ^^^^ I I logl^-^/rVi(a^)^fc(y) + X]'^^' / + 
,=1 k=i+i i=i 



c := lim \ \ogZr, 



Corollary 3.2. Whenever the joint density of the eigenvalues is as in i \3.2\\ and the rate func- 
tion I has a unique minimis er , we obtain under (al) and (a2) a strong law of large numbers, 

P(L„ ^ /i*) = 1. 

Remark 3.3. Hence the vector of empirical measures of a random points from the Angelesco 
ensembles tends to the vector of nonrandom measures almost surely weakly. This solves the 
question posed in [111 Section 5.2]. Remark that recently, in a strong law for Angelesco 
ensembles was established, applying the notion of Fekete points as well as the Bernstein-Markov 
inequality. We obtain a full LDP. We would also be able to consider Nikishin ensembles with 
p > 2 weights, see [H] and references therein. This is because the determinantal structure 
of the joint density of the eigenvalues [HI (4.14)] consists of Vandermonde-like products and 
hence the techniques of our proof of Theorem 12.11 can be adapted. Nikishin interaction arises 
in the asymptotic analysis of eigenvalues of banded Toeplitz matrices as well as in a two-matrix 
model, see [H] Section 5.4]. 



4. Proofs 

Now we present the proof of Theorem 12.11 It is a generalisation of arguments used for the 
proofs of [H Theorem 2.6.1] and [HI Theorem 4.1]. The proof of the upper bound is quite 
similar to the proofs in the latter references. In the proof of the lower bound, a coarse graining 
argument is more involved. To overcome the singularity of the logarithm at certain points is 
the most delicate part. 

Let us define F : S x S ^ 1 by 

F{x,y) := -— (log|x-y| + log -/|) - -{\ogw{x) + hgw{y)), 
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where we set F{x, y) = oo if = or if {x, y} fl Af{w) ^ 0. Note that due to our definition 
of E, x^ = y^ corresponds to x = y. Let denote the truncated version of F{x,y), for 

M > 0: F^^{x^y) := F{x,y) A M. Furthermore, we define the functions F„ : S x S — )■ M, 

Fn{x,y) := -- [—^j (log k - Z/l + log l^^^ " /|) " (log (a;) + logw„(?/)), 

again with Fn{x, y) = oo if x^ = y^ or if {x, y} fl M{wn) 7^ and their truncated versions 

F^'{x,y) ■.= Fn{x,y)AM, M > 0. 

Observing that 

p(n) ^ ^ p(n) 

n V'logw„(xi) = — - (logti?„(xi) + logty„(a;j)) + — - logti?„(xi), 

we can deduce from (12.11) and the definition of Fn the following identity, where we abbreviate 

lj:p(n){xi, ..,Xp(^n)) by l£p(„)(x), 

1 f 2n2 n ^^"^ 1 

g„(Xi, . . . = — exp<J ——y^Fn{Xi,Xj) + ——y^logWn{Xi) > l2P(n)(x). (4.1) 

1 P(n)2 ^ p(n) ^ J 

One key ingredient for the proof of the upper bound will be the following lemma, which also 
provides that the rate function is well defined. 

Lemma 4.1. 

(i) For any M > 0, F*^(x, y) converges to F*^(x,?/) uniformly as n 00. 

(ii) F is bounded from below. 

Proof. Since log \x — y\ < log(|x| + 1) + log(|?/| + 1) holds for any x,y eM., it implies 



Fn{x,y)> ^^"^ 



( pin) \ / p(n) 

[i\x\ + i)i\x'\ + i)] " t^„(x)j +iog( [(|y| + i)(|/| + i)] " My) 

We will show that F*^ = F^^ = M on some specified sets and then deal with the complement 

( p(n) \ 
[(|x^| + l)(|a;| + 1)] Wn{x) j is bounded 

from above: In case of x G [—1, 1] fl S it turns out that 

log (^[(|X^| + 1)(|X| + 1)] Wn{x)^ < log (4^^) + \og{Wn{x)). 

For M > we can take Ui as large, so that (al.l) is applicable and — n\ < e and 
\wn{x) —w{x)\ < e V n > rii and V a; G [—1, 1]. For each v G M{w) a 5^^\i > and V a; G [—1, 1] 
with |x — z/| < Sl^lj, we now find log (^4^^ + log(w„(x)) < log(4''+^) + log(w(x)) < —M. 
Whereas for |x| > 1, we observe that for n > ni as above the following holds 



log (|a;'^| + l)(|a;| + 1) " w^ix) 



p(") 



) ) < log U^+^\x\(o+^^(^+^^ sup Wm{x)j . (4.2) 

/ \ m>n / 
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Assumption (a2) provides that for each M > there exists n2 E N and Rm > such that for 
> Rm and all n > ^2 we have that 



log U-+^\x\'^<^+^^^^+^) sup 

\ m>n 



w^x) < -M. 



In case of < Rm reasoning as in the first part yields for each M > and u G N'lw) the exis- 
tence of > 0, such that for x with |x — z/| < 5^,^],^ we also get log |a;|(^+^)^ Wn{x)j < 

—M. Thus, for each M > there exists uq := max{r;,i, 77-2}, Sm,u '■= ''^^''^{^^I,u^^^M,u} > 
Rm > such that F„(x, y) > M holds for all n > no on 

AM:={\x\y\y\> Rm}^^ U {\x-u\A\y-u\<5,,M}- (4.3) 

From (14.31) we see that A\,j- is compact. Moreover, we find a constant Cm > 0, depending on M, 
such that Fn > M and F > M W x,y E Bm '■= {{x, y) G : |x — y| < Ca/}, and therefore, a 
constant < Cm < 00 exists with max^^yg^c^j— log \x — y\, — log — y^l} < Cm- This implies 
the existence of yet another constant < Ci^m < 00 such that 

max {|log|x -yll, |log|x^ -y'^ll} < Ci^m- 

x,yGBl,nAli 

Remember that x^ = y^ only holds if x = y. Furthermore, we introduce the notation || 
/ II sup^g£) I /(a;) I, which restricts the supremum norm to a set The continuity of Wn 
and the compactness of Alj yields that || logw„(x) + logWn{y) ||oo"< C2,m, for a constant 
< C2,M < 00. Thus, for any given 77 > we find on D := A%i \ Bm an G N, such that 
V n > n 

II F.^{x,y)-F{x,y) \\Z<V- 

Hence, we have established the uniform convergence of F^^ to on A%i \ Bm-, whereas on 
Am U Bm we have F^^ = M. That F*^ = M holds on Am U Bm can be shown along the same 
lines and (i) is proven. Besides F > M on Am ^ Bm, we also know that F is real- valued and 
continuous on the compact set A1^ D B%j and therefore F is bounded from below, which yields 
(ii). □ 

Before proving the LDP for [Qn o L^^)„, we consider first the finite measure 

where Z„ is the partition function of the Lebesgue density belonging to Q„, see (12. ip . and derive 
a LDP for (P„oL-i). We set 

H[ii) ■= j Fd^i^^ H^\^) := j F^'d^i'^^ 

and obtain well defined maps on A^i(S) (due to Lemma bounded). We claim the following: 

Lemma 4.2. H is a good rate function that governs the LDP for {Pn o L~^)n. 
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Proof. Exactly as the proof of [SI Le. 3.5], since Lemma [4.11 assures that F has the necessary 
properties to mimic the proof of [SI Le. 3.5]. □ 

Moving on, we will show that (P„ o L~^)n fulfils a weak LDP. 

4.1 Proof of the upper bound The main obstacle will be to overcome the singularities of 
the function F. We start by observing that for A := {(a;, G : x = y}, 

Ln ® Ln{^) = (4.4) 
p[n) 

holds mj5p(n) -almost surely, because of the a.s. distinct eigenvalues under the product Lebesgue 
measure mj^p{n). From (14. 4 p it follows that Q„-a.s. 

J J^^ F^'{x,y)Ln{dx)Ln{dy) = J J F^' {x,y)Ln{dx)L^{dy) - (4.5) 

For a Borel set A in A^i(S) let us define A as follows A := [x e S^^") : L„(x) G A}. Due to the 
symmetry of F„ in its arguments we find 2 J2i<j Fn{,Xi, xj) = X^i^^j Fn{,Xi, Xj), and since mj^ptn) 
is the product Lebesgue measure on R^'^"^ we obtain via Holder's inequality 

P„(L„ eA) < fy exp |^logw„(t)| m^idtU 



1/2 



= (I)x(II). 

First, we show lim^^oo log (I) = 0. Since Wn is continuous, hence bounded on compact sets, 
the key observation is that (a2) implies 



/- 



2n 



w;„(t)j'(") mM(rft) = / Wn(t)H") w . mK(rft) 



for suitable K,ri > and n large enough. Concerning (II), we find that for any M > 

1/2 



(II) ^ {/^exp(-2„=(L„(.r(F„"')-^) 



m]gp{n) (dx) 



< <'exp( -2nM inf /i^^^Ff)- 



M 



MSA y 



1/2 



exp ( -n^ inf /x®'(Ff )^ exp f ^ ) . (4.6) 
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The first part of Lemma 14.11 yields, 

lim ( inf )^ = inf fi'^^F^'). 

We have thus shown that for any Borel set A C A^i(S) one has 

limsup^logP„(L„ e A) < - inf / [ F^\x,y) fi{dx)n{dy). (4.7) 



The inequality (14. 7p allows us to proof the following lemma. 
Lemma 4.3. (P„ o L~^)n is exponentially tight. 

Proof. Very similar to ^ P-12] so that we omit the proof here. □ 
Now we choose for A the set -B(/i, 6) = {u E A^i(S) : c/(z/, /i) < 26} with 6 > and 



fd/j,- / /rfi/ 



d{fi, p) = sup 

where the supremum is taken over all Lipschitz functions / for which the sum of the Lipschitz 
constant If and of the uniform bound || / ||oo is less than or equal to 1. That distance is 
compatible with the weak topology, see P p. 356]. Since /x h-)- H^\ij) is weakly continuous, 
from (14. 7p we obtain for any ji E A^i(S), 

1 E3 f f 

inf lim sup — log P„(L„ G -B(/i, 5)) < —inf inf / / F^Ux,y) u(dx)u(dy) 

Finally, letting M go to infinity, we obtain the following upper bound, 

inf limsup^logP„(L„ E B{fi,6)) < -H{fi). (4.8) 



4.2 Proof of the lower bound Turning to the lower bound, we take P(/i, 6) as above and 
for every E A^i(S) we will show a lower bound, 

inf liminf 4logP„(L„ E B{fi,6)) > -H{fi). (4.9) 

S>0 n-i>oo 

As in [2] and [8] we can assume w.l.o.g. that 

(i) fi has no atoms, 
(ii) S := supp(/i) is a compact subset of S with S fl {Af{w) U {0}) = 0. 

The main idea of the proof of the lower bound, is to localise the eigenvalues in small sets, and 
benefit from the speed n^, which provides that the small volumes of these sets can be neglected. 
We will divide S as follows. For j = 1, . . . ,p{n) let = ^j"-' the ^^""p^^'^ quantile of /i and 
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set = (^p(n), . . . ,^i), ^p(n)+i := inf 5 and ^0 = ^1 + 1- Therefore, ^(] - oo,^p(„)+i]) = 
oo[). Due to assumption (ii), we have 

-00 < ^p(„)+l < ^p(n) < • ■ ■ < ^1 < ^0 < OO- 

For 5 > 0, t G write 

• 7r„(t) := {i = 1, . . .,p{n) : U > 0}, z/„(t) := {1, . . .,pin)} \ 7r„(t), 
. 4(5):={z = 1,...,pH: |d"^-d?il<5}, 

• X]"^(t, (5) := [t, - S, tj + 5] n S, J = 1, . . . ,p(n), 



+ forjG7r„(t), 
[tj-5,tj] for jeun{t), 



We may assume that 6 is fix with < 6 < Set 

ifi'f^ := v?^"'''^ := inf{w„(x) : x G [^^ - 6, + S]U [^j+u ^j-i]}, j = 1, • • • ,p( 



n] 



and analogously define ipj when w;n(a^) is replaced by w{x) in the above definition. Write ip„ 

4' 



resp. ip for the step function which equals (/^^"^ resp. on and is zero elsewhere 



p{n) p{n) 

i=i j=i 



Moreover, we have: 



Lemma 4.4. 3 G N s«c/i ^/iat V n > no : C {x G S^^") : L„(x) G 5(/i,2(5)}, 

w/iere B{n,26) = {u e Ali(S) : d(z/,/i) < 25} with 



d(fi, v) = sup 

{fLipschitz: i/ + ||/||oo<l} 



fdfi- / /dz/ 



as above. 



Proof. Let y = {yp{n), y,{n)) G S). Thus, G [e. - <5, 6 + <5] and 

. pin) 

d{LM, L„(e("))) < sup \fiy^) - < \y^ - 61 < ^• 

i/ + ll/l|oo<l «=1 

Applying [21, Le. 3.3] finally yields /i) < + /x) < 25. □ 
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Now for n > uq, we get from Lemma HiH and I„(^, S) D J„(^, S) the inequality 



p„(L„ei?(/.,25)) > z„Q„(i„(e,5))>nM"0 / Hk'-s 



Xi - Xjl mp(„)(da;). 

(4.11) 



Introducing the notation 



we focus on the last integral, 



{x e M'^lxi > X2 > . . . > Xd} 



> 



/ TT |x- - x^l |xi - Xjl mp(„)((ix) 

Y\\ixi + ^i) - (xj + 0)1 \ixi + ^i)' - (xj + 0)'| mp(„)(da;) 

-5,0]#''"(«x[0,<5]#''"(«) 
([-<5,0]#''"(«)x[0,5]#''"(«))nIRP(")'+ 



([-(5,0]#''"(«)x[0,<5]#'^"(«))nRP( 



|xj Xj + 01 I (Xj + ^j) + 



.i<i-l 



p(n)-l 



|Xj — Xj+i + — 0+1 1 I + 0)^ ~ (3^1+1 + ^i+lY 



i=l 



m 



„)((ix). 



On ]RP(")'+, the following inequalities hold, 

Y[ \xi-xj+^i-^j\> Y[ 16-01 

and as it can easily be seen (e.g. from ab < {a + 6)^, for a,b > 0) 

p(n) — l p{n)~l 



(4.12) 



(4.13) 



JJ |xi -Xi+i +0 -O+il > n -a;i+i|2|0 -O+il^' 



j=l 



Next, we will show that 



(xi + 0)' - {xj + 0)1 > + - a;,' - O'l > 



(4.14) 



(4.15) 



where the modulus can be omitted since both terms are greater than on Mp("-)'+ due to our 
choice of the 0- If > > ^ ^^"^ Xi > xj > 0, we find that 



Xi~\~^i) (•^i~l~0) •^i'^^i •^i 



>o 



e-k _ kfd-k 
Xj 



In case of > > ^^"^ > Xj > Xj, we first recall that this only can happen while is odd. 
Nevertheless, (I4.16P also holds: Within the sum, either ^ — is odd and k is even or vice versa. 
Hence, if k is even > i'l'^x^i > O^-'^xf > O^"''^, while the other case follows in an analogous 
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way. The last case to consider is > > and Xi > > Xj. This case also only comes up for 
9 being odd, and we find xf^f"'^ > and > since either k is even or 9 — k. Thus, we 

have proven fl4.15p and this can on one hand further be used to obtain either 

x^+^!-x^-^>ix^-x^rie!-^r (4-17) 

or on the other hand yield 

x^ + ef-4-e?>ef-^, (4.18) 
since xf > on Mp(")'+ for 6* G N. 

Putting together IKT^ . KIM . iKTf} and KIM , we can bound fl412|) from below by 

p(n)-l ^ 
i<j—l i=l 

/ n ~ ^i+l)^ i^i ~ ^i+l)^] ^Pin){dx). (4.19) 

^([-(5,0]#''"(«)x[0,<5]#'^"(O)nMP(").+ L J 

Applying the mean value theorem for 9 > 2 (while for 9=1 the following inequality will be an 
equality) now yields x^ — x^^i = {xi — Xi+i)9(^^^, with ( E [xi+i,Xi]. Since \xi\ < 6 W i and 
6 < , we easily see that < 9(^^^ < 1. Therefore, < xf — xf , ^ < Xj — Xj+i and we obtain 



> 



([-<5,0]#''"«)x[0,<5]#'^"(«))nRp("). ■ j^]^ 

p(n) — 1 

n K -^m)"^p(")(^^) 

([-<5,0]#''"(«)x[0,<5]#'^"{«))nRP(")'+ j^;^ 

p(n) — 1 



/ n {xi ~ Xi+i)\\\xi\' mp(n){dx) (4.20) 



i=l 1=1 



> / n (xi-Xi+i)mp^^){dx). (4.21) 

^ J([-<59,0]#''"(Ox[0,5e]#'^n(f))nRP(") ' 



The equality fl4.20p is a simple application of the transformation formula, while the inequality 
( 14.2 ip holds, because 9 is an integer and 5 < 1, so that |xj|^ > 1. For the last integral, we 
substitute = Xp(„), = Xj_i — Xi,i = p{n), . . . , 2 and obtain 

p pH-I p pin) 

/ IT {xi - Xi+i) mp(^n){dx) = / ]TMimp(„)((ix), 

^([-<5«,0]#''"(«)x[0,<59]#'^"(«))nRP(")'+ ~~l -^^^1=2 

P{n) /^\p{n)-l / ^0 ^ 2{p(n)-l)+l 
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where [/„ = |m = (Mp(„), . . . ,Mi)|Mp(„) e [-5^, 5^], G [0,(5^],i = 1, . . . ,]?(n) - 1 : | Ya=i < '^^} 
From fHTT]) . fICTD and fH:^ we get, 



^ log P.(L„ei?(/., 25)) > ^ Elog([te-e.)(C^-e,^)]) (4.23) 

+ (^) ^ E '"S ([(«.-«.«)(«.' (4.24) 
+ -El^g^S"^ (4-25) 



i=l 

;^log (2(^)'^"V^^^^(log5^-logMn)). (4.26) 



It is easily seen that f l4.26p converges to as n tends to infinity. 
Now, for (i23D and (i21D observe that for 6^ = 1, resp. 6* e N, 



log(y^ — x^)fi{dx)^{dy) — / log \y^ — x^\fi{dx)fi{dy) 



^ / log(?/^ - x'^)fi{dx)fi{dy) 

+- V / log 1/ - sV(f^a;)/i(dy) 



=1 >^{a;,y)eK,+i,5«]' 



Because was defined to be the infimum of the support S and we look at measures without 

atoms, we observe that ^p{n) — ^ ^p(n)+i = inf 5 for n — )■ oo. Thus, we find 

lim / log \y^ — x^\^{dx)fi{dy) = / log \y^ — x^\fi{dx)fi{dy) = 0. 



We have shown that the limit of (14.231) and (14.241) can be bounded from below in the following 
way, 

j j log(l/^ - x^)^i{dx)n{dy) = \j^ log(l/^ - x^)^i{dx)^i{dy) 

s \ f^;,^ E - ?|«) + ^ 'E - «^..)) ■ (4-27) 
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As to fl4.25p . observe that 

p{n) 



1 v^i (n) pin) f 
n ^—^ ■' n J 



n 

3 = 1 



where ipn was defined in fl4.10p . Denote by / a Lipschitz constant of \ogw on S according to 
(al.3). For > write 

:= max{| \ogw{x) — logw(y)| : \x — y\ < rj}. 

Note that < /?7, again because of the Lipschitz continuity (al.3) and since N{w)r\S = 0. Next, 
define M{n,6) := Ujg/„(5) [^j+ii ^j"^] ^ ■= max{logw(x) : x G 5}— min{logw(x) : x G 5}, 
where assumption (ii) provides that C is finite. Let e > 0. Since for all n and j > 1 one has 
^Y^]^ = l/p(^)) since 6 is fixed, one has fi{M{n, 6y) < e for large n. Now let n be 
large enough such that one also has || logw„ — logt(;||oo < e, which is possible because of (al.2). 
Then 

/ I log '?/'„, — log w I (i/i < / {| log'?/' — log w| + e} (i/i 

+ ^ / {| log'?/' - logw| + e} (i/i + e 
< p(n)-^ (/2^ + e) + (C + e)e < 2W + (C + e + l)e 



This yields 



lim sup 



— log yj^"^ — log u; d/i 



0(5). (4.28) 



Putting together fl4.23l) - fl4.28p yields the lower bound ( 14. 9p . whenever the empirical measure 
lies in a ball around some measure n in A^i(S). 

4.3 Extending the weak LDP into a full one As S is a closed subset of the Polish space 
M, A^i(S) is a Polish space, cf. [15, Thm. 6.2] and since the B{fi,6) are a base for the weak 
topology on A^i(S), we can apply [6l Thm. 4.1.11], which provides a weak LDP for (P„oL~^)„ 
on A^i(S) with rate function H and speed n^. Because (P„ o L~^)n is also exponentially tight. 
Lemma 1.2.18 in [6] extends the weak LDP into a full one. But our original aim was to derive 
a LDP for (Qn o L~^)n- Setting A = G = A^i(S) in the lower and upper bound at the LDP 
for (Pn o L~^)n, we obtain 

lim 4log^n = - inf [ Fdfx'^^ 

n^oo /ieA4i(S) J 

The right hand side has been shown to be < +oo by Lemma UT](ii) and establishes (12. 3p . Now, 

^logQ„(L, eA) = \ (log PniLn G A) - logZ,) 
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for any Borel set A in A^i(S). Hence Theorem 12.11 is proven. 

5. Sketch of the proof of Theorem LH.ll 

The Vandermonde-hke products A(X) and A(X, Y) in (13. 2 p lead to the observation, that to 
overcome the singularity of the logarithm at several places is the most delicate part. Therefore, 
we will skip most of the technical details of the proof worrying about too many replications of 
our arguments. Instead we give a sketch to be able to discover the rate function. Consider a 
neighborhood O of fi = {fii, ... , /ip) G Mi{Tp^fi). With O := {x G M" : -L„ G 0} we obtain 

Qn{Ln e O) = I ■ ■ ■ I exp{F{x))dxi . . . dxn 



Zn J JO 

with 

P rii P ni 

n^) = E E - + E E E - ^i') - ^ E E ^^(4^^ (s.i) 

i=l l<j<k<ni ^<i<j<P k=l 1=1 i=l k=l 

Very roughly speaking, we consider the case when O goes to a point ^. Then for x G O the 
first summand in (15.11) can be approximated by 

E"^ / hog\x -y\'^diJ.,{x)dij,i{y), 

the second summand by 



^ UiUj \og\x - y\diJ,i{x)diJ,j{y) 

l<i<j<p 

and the third summand by n Yl'i=i J yi{x)dfti{x). With (13. 3p one gets 
^logg„(L„ G O) 

P ^2 



~ E^ / hog\x - y\^dfii{x)dfii{y) + ^ r.rj log \x - y\dfii{x)dfij{y) 

i=l l<i<j<P 

+ E'^^y Vi{x)d^ii{x) + — log Zn. 
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